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ABSTRACT 

Simulations have indicated that most of the escaped Lyman continuum photons 
escape through a minority of solid angles with near complete transparency, with the 
remaining majority of the solid angles largely opaque, resulting in a very broad and 
skewed probability distribution function (PDF) of the escape fraction when viewed at 
different angles. Thus, the escape fraction of Lyman continuum photons of a galaxy 
observed along a line of sight merely represents the properties of the interstellar medium 
along that line of sight, which may be an ill-representation of true escape fraction of the 
galaxy averaged over its full sky. Here we study how Lyman continuum photons escape 
from galaxies at z = 4 — 6, utilizing high-resolution large-scale cosmological radiation- 
hydrodynamic simulations. We compute the PDF of the mean escape fraction ((/ esc ,id)) 
averaged over mock observational samples, as a function of the sample size, compared 
to the true mean (had you an infinite sample size). We find that, when the sample size 
is small, the apparent mean skews to the low end. For example, for a true mean of 6.7%, 
an observational sample of (2,10,50) galaxies at z = 4 would have have 2.5% probability 
of obtaining the sample mean lower than (/ eS c,iD) =(0.007%, 1.8%, 4.1%) and 2.5% 
probability of obtaining the sample mean being greater than (43%, 18%, 11%). Our 
simulations suggest that at least ~ 100 galaxies should be stacked in order to constrain 
the true escape fraction within 20% uncertainty. 


1. Introduction 


A fraction of the Lyman continuum (LyC) photons generated by young massive stars is believed 
to escape from the host galaxies to enter the intergalactic space. This is a fundamental quantity to 
determine the epoch and pace of cosmological reionization, provided that the universe is reionized 
by stars (e.g., Gnedin 2000 Cen 2003). After the completion of cosmological reionization, it plays 
another important role in determining the ultra-violet (UV) radiation background (on both sides 
of the Lyman limit) in conjunction with another major source of UV photons - quasars - that 


progressively gains importance at lower redshift (e.g., Faucher-Giguere et al. 2008 Fontanot et al. 


2014). 


Observations of star-forming galaxies at high redshifts (z ~ 3) suggest a wide range of the 
escape fraction of ionizing photons. While only a small fraction of LyC photons (< a few percent) 
escapes from their host galaxies in the majority of the Lyman break galaxy samples, a non-negligible 
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number of them (~ 10%) shows high levels of LyC flux corresponding to (/ eS c,iD) ~ 10% (£ 

ha p ley 

et al. 

2006 

Iwata et al. 2009; 

Nestor et al. 

2011 

2013 

Mostardi et al. 2013). Cooke et a. 

(2014) 


claim that the mean escape fraction may be even higher ((f esc , ID) ~ 16%) if the observational 
sample is not biased toward the galaxies with a strong Lyman limit break. It is not well understood 
quantitatively, however, what the probability distribution function (PDF) of the LyC escape fraction 
is and how a limited observational sample size with individually measured escape fractions can be 
properly interpreted, because of both possible large variations of the escape fraction from sightline 
to sightline for a given galaxy and possible large variations from galaxy to galaxy. The purpose 
of this Letter is to quantify how LyC photons escape, in order to provide a useful framework for 
interpreting and understanding the true photon escape fraction given limited observational sample 
sizes. 


2. Simulations 


To investigate how LyC photons escape from their host halos, we make use of the cosmological 
radiation hydrodynamic simulation performed using the Eulerian adaptive mesh refinement code, 


(Teyssier 2002 

Rosdahl et al. 2013 

ver. 3.07). The reader is referred to 

Kimm & Cen 

(2014 


, the FRU run) for details, where a detailed prescription for a new, greatly improved treatment of 
stellar feedback in the form of supernova explosion is given. Specifically, the new feedback model 
follows the dynamics of the explosion blast waves that capture the solution for all phases (from early 
free expansion to late snowplow), independent of simulation resolution and allow for anisotropic 
propagation. 


The initial condition for the simulation is generated using the MUSIC software (Hahn & Abel 


2011), with the WMAP7 parameters (Komatsu et al. 2011): (h2 m , Da, 1%, h, as, n s = 

0.272,0.728,0.045, 0.702, 0.82, 0.96). We adopt a large volume of (25Mpc/h) 3 (comoving) to include 
the effect of large-scale tidal fields on the galaxy assembly. The entire box is covered with 256 3 
root grids, and high-resolution dark matter particles of mass M c \ m = 1.6 x 10 5 M 0 are employed in 
the zoomed-in region of 3.8 x 4.8 x 9.6 Mpc 3 . We allow for 12 more levels of grid refinement based 
on the density and mass enclosed within a cell in the zoomed-in region to have a maximum spatial 
resolution of 4.2 pc (physical). Star formation is modeled by creating normal and runaway particles 
in a dense cell (rin > 100cm~ 3 ) with the convergent flow condition (Kimm & Cen 2014, the FRU 
run). The minimum mass of a normal (runaway) star particle is 34.2M 0 (14.6M 0 ). We use the 
mean frequency of Type II supernova explosions of 0.02 M 0 _1 , assuming the Chabrier initial mass 
function. Dark matter halos are identified using the HaloMaker (Tweed et al. 2009). 


Eight consecutive snapshots are analyzed at each redshift (3.96 < z < 4.00, 4.92 < z < 5.12, 
and 5.91 < z < 6.00) to increase the sample size in our calculations. At each snapshot there are ~ 
142, 137, and 104 halos in the halo mass range of 10 9 < M v - U < 10 10 M 0 , and 15, 10, and 7 halos 
with mass M vir > 10 10 M 0 . The most massive galaxy at z = 4 (5, 6) has stellar mass of 1.6 x 10 9 M 0 
(6.0 x 10 8 , 2.5 x 1O 7 M 0 ), and host halo mass 8.8 x 1O 1O M 0 (5.2 x 10 10 , 4.1 x 10 11 M 0 ) 


The escape fraction is computed as follows. We cast 768 rays per star particle and follow their 
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propagation through the galaxy. Each ray carries the spectral energy distribution (SED), including 
its LyC emission, determined using STURBURST99 (Leitherer et al. 1999), given the age, metallicity, 


and mass of the star particle. The LyC photons are attenuated by neutral hydrogen (Osterbrock 


& Ferland 2006) and SMC-type dust (Draine et al. 2007) in the process of propagation. For a 


conservative estimate, we assume the dust-to-metal ratio of 0.4. We also simply assume that dust 
is destroyed in hot gas (T > 10 6 K). We note that attenuation due to dust is only significant in the 
most massive galaxy (M star = 1-1 x 10 9 M o , r<j = 0.58) in our sample. The second most massive 
galaxy (M star = 3.6 x 10 8 M e ) shows Td = 0.29, meaning that it reduces the number of photons by 
only < 30%. Given that the dust-to-metal ratio is even smaller than 0.4 in low-metallicity systems 


(Lisenfeld & Ferrara 1998; Engelbracht et al. 2008 Galametz et al.l 2011 Fisher et al. 2013), it is 


likely that the attenuation by dust is even less significant in our simulated galaxies. We define the 
true escape fraction of the galaxy as the ratio of the sum of all outward fluxes at the virial sphere to 
the sum of the initially emitted fluxes of all stellar particles in the galaxy; we shall call this / eS c,3D- 
In addition, an observer at infinity at a random point in the sky of the galaxy collects all LyC fluxes 
and defines the escape fraction along that particular line of sight; this is called / esc ,iD- 


3. Probability Distribution Functions of LyC Photon Escape Fraction 


ft is useful to give a qualitative visual illustration of how LyC photons may escape from galaxies 
at z = 4. Figure [T| shows three examples of an all-sky map - the sky an observer sitting at the center 
of the galaxy would see - of the neutral hydrogen column density. We note that 8 dex of dynamical 
range is plotted and recall that at the Lyman limit a neutral hydrogen column density of r^j 10 17 cm 2 
would provide an optical depth of ~ 1. As a result, LyC photons can only escape through highly 
ionized or evacuated “holes” indicated by dark blue colors on the maps and the transition from near 
transparency to very opaque is fast. This indicates that the escaping LyC photons are dominated by 
those that escape through completely unobscured channels and the amount of escaped LyC photons 
for a given galaxy depends strongly on the direction. Moreover, it is evident that, in addition to 
large variations from position to position on the sky for a given galaxy, there are large variations of 
the overall column density structures from galaxy to galaxy. For example, the galaxy in the top-left 
panel shows no transparent sky patches at all, which is typical for galaxies during times of intense 
starburst as shown in Kimm & Cen (Figure 4 2014). On the other hand, the galaxy in the bottom 
panel has large swaths of connected transparent patches that cover nearly one half of the sky, typical 


for galaxies at periods following the blowout of gas subsequent to intense starburst (Figure 4 Kimm 


& Cen 2014). This qualitative behavior is also found earlier in independent simulations by Wise & 


& Cen 

2014) 

Cen 

(2009 

)• 


Let us now turn to more quantitative results. Figure [2] shows the probability distribution of 
the apparent escape fraction for massive halos (top) and less massive halos (bottom) at z = 4 (left 
column) and z = 6 (right column). Black histograms show the distribution of the true (3D) escape 
fraction of each sample (i.e., from the viewpoint of the overall intergalactic medium), while red 
histograms show the PDF of the apparent escape fraction (i.e., from the point of view of observers 
placed at a far distance). Note that the distribution of the true escape fraction is noisier than that of 
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Fig. 1.— shows three examples of all-sky maps - the sky an observer sitting at the center of the 
galaxy would see - of the neutral hydrogen for most massive (M vir = 7.8 x 1O 1O M 0 , top left panel), 
second massive (6.1 x 10 10 M o , top right panel), and a smaller halo (1.8 x 1O 9 M 0 , bottom). The 
observer is placed at the center of the halo. Note that the actual escape fraction presented later 
is computed by ray-tracing LyC photons of all stellar particles spatially distributed through the 
clumpy interstellar medium until escaping through the virial sphere. The true escape fraction of 
LyC photons of these halos are 5.4%, 12%, and 5.0%, respectively. 

the apparent escape fraction due to the smaller sample size for the former, because for (3D) escape 
fraction each galaxy is counted once but for the apparent escape fraction each galaxy is sampled 
many times. In terms of the mean escape fraction, there is a trend that, at a given redshift, the 
galaxies embedded in more massive halos tend to have a lower mean escape fraction. 

There is also a weak trend that the escape fraction increases with redshift. For example, the 
true (3D) median escape fractions are (7.0%, 9.5%) for the halos of masses (> 10 10 ,10 9 — 10 10 ) M 0 , 
respectively, at z = 4; the true (3D) median escape fractions are (8.8%, 29%) for the halos of 
masses (> 10 10 ,10 9 — 10 10 ) M 0 , respectively, at z — 6. Upon a close examination we suggest that 
the redshift dependence can be attributed, in part, to the following findings. At a given halo mass, 
the specific star formation rate decreases with decreasing redshifts at 4 < z < 6. As star formation 
becomes less episodic at lower redshifts, it takes longer to blow out the star-forming clouds via 
SNe. Consequently, a larger fraction of LyC photons is absorbed by their birth clouds. We also 
find that the specific star formation rate does not change notably at z > 6 while the mean density 
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Fig. 2.— shows the probability distribution of the apparent escape fraction for massive halos (top 
panel) and less massive halos (bottom panel) for z = 4 (left column) and z = 6 (right column). Black 
histograms show the distribution of the true escape fraction of each sample, while red histograms 
show the PDF of the apparent escape fraction. The median of the distribution is shown as arrows. 


of the halo increases with redshifts, explaining an opposite trend found in Kinim & Cen (2014) 


at 7 < z < 11. The predicted median escape fraction for halos with 10 9 < M halo < lO n M 0 at 


4 < z < 6 is generally smaller than the previous studies (10 — 30% Razoumov & Sommer-Larsen 


2010 Yajirna et al. 2014), although the trend at z < 6 is broadly consistent with Razoumov & 


Sommer-Larsen (2010). Nevertheless, it is prudent to bear in mind that the sometimes conflicting 


results and trends among studies with respect to redshift may be in part due to still limited galaxy 
sample sizes. 

Figure [3] is similar to Figure [2j except the galaxy sample is subdivided according to their star 
formation rates (SFRs). We note that this division according to SFRs introduces subtle degenera¬ 
cies. For example, a lower SFR does not necessarily correspond to a less massive galaxy; instead, a 
lower SFR may correspond to the phase of a galaxy between two star formation bursts. The signif¬ 
icantly higher escape fraction for the lowest SFR bin (bottom panels) is, to the most part, due to a 
post-starburst phase when the interstellar medium has been cleared out by the preceding burst and 


SFR has abated, as noted in Kimm & Cen (2014). Thus, if we do not consider the lowest SFR bin, it 
seems that the mean escape fraction does not strongly depend on SFR at z = 4 — 6. The escape frac¬ 
tion from the most actively star-forming galaxy sample with 0.3 < SFR < 10M e /yr be compared 
with that of Lyman alpha emitters or faint LBGs. Our simulations suggest that the median / eS c,iD 
of the sample is 5.1%, which is consistent with / es c ,id = 5 — 15% inferred from narrow band filter 
imaging observations of 91 LAEs (Mostardi et al. 2013). We note that, given the wide distribution 
of the simulated apparent escape fraction (0.1% < / esc ,iD < 23% or 0.001% < /esc, id < 47% for the 
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1 and 1.5 a range, respectively), onr results are also compatible with the individual detection of 
LyC fluxes from 7 LBGs (Iwata et ah 2009, 5.5 < / esc , id ^ 55% for the intrinsic UV to LyC flux 
ratio of 3.). 



Fig. 3.— is similar to Figure[2j except the galaxy sample is subdivided according to their star forma¬ 
tion rates, SFR=0.3 to 10M o /yr (top panel), 0.01 to O.3M 0 /yr (middle panel), and < 0.01M Q /yr 
(bottom panel). 

Evidently, the distribution of the apparent LyC escape fraction is very broad and skewed toward 
the lower end. The reason for this behavior is understandable. In the case of the galaxies with low 
fesc,3d values, the LyC photons escape normally through transparent holes with small solid angles. 
Since not all of these holes are seen to an observer, the distribution of f esc ,id for individual galaxies 
tends to get skewed toward the lower end of the distribution. As a result, the medians of the two 
distributions, shown as arrows in Figure [2j are about a factor of ~ 2 smaller than the mean. More 
importantly, it suggests that an observational sample of limited size may underestimate the true 
mean escape fraction. The top two panels of Figure [4] show the probability distribution function 
of the apparent mean for a given observational sample size Attack for the high mass (top) and low 
mass (bottom) sample, respectively. We compute the apparent mean of a sample of galaxies using 
LyC photon (or SFR)-weighted mean escape fraction, which is exactly equivalent to stacking the 
galaxies. The bottom two panels of Figure [4] are similar to top two panels in Figure [4], for the 
subsamples with different star formation rates. What we see in these figures is that the probability 
distribution is rather broad. It is thus clear that it is not a robust exercise to try to infer the 
mean escape fraction based on a small sample (< 10) of galaxies, whether individually measured or 
through stacking. 

Table 1 provides a quantitative assessment of the uncertainties, which shows the 1 and 2 a 
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probability intervals of fractional lower and upper deviations from the true mean escape fraction. 
Some relatively mild trends are seen that are consistent with earlier observations of the figures. 
Specifically, the convergence to the true mean escape fraction in terms of sample sizes is faster 
towards high redshift, towards higher halo mass, and towards higher star formation rates. Let 
us take a few numerical examples. We see that with a sample of 50 galaxies of halo mass in the 
range of (10 10 — 10 11 ) M Q at z = 4 the 2cr fractional range of the escape fraction is 58% to 159%, 
which improves to a range of 68% to 140% when a sample of 100 galaxies is used. Note that the 


observations of Mostardi et al. (2013) have 49 Lyman break galaxies and 91 Lyman alpha emitters 
at z ~ 2.85. At z = 6 for the (10 10 — 10 11 ) M 0 halo mass range, we see that with a sample of 
20 galaxies, the 2cr fractional range of the escape fraction is 59% to 161%, comparable to that of 
a sample of 50 galaxies at z = 4, as a result of benefiting from the faster convergence at higher 
redshift. On the other hand, at z = 5 for the (0.3 — 10) M 0 yr _1 star formation rate range, the 
2a fractional range of the escape fraction is 56% to 163% with a sample of 20 galaxies, which is 
improved to 71% to 137% with a sample of 50 galaxies. 

Finally, we note that the actual observed Lyman continuum escape fraction has additionally 
suffered from possible absorbers in the intergalactic medium, primarily Lyman limit systems. Since 
the background galaxy and the foreground absorbers are physically unrelated, we may consider the 
effects from the internal factors in galaxies and those from the intergalactic medium completely 
independent. Thus, in this case, assuming no knowledge of the foreground absorbers, the overall 
distribution would be the convolution of the two, resulting a still broader overall distribution than 
derived above considering internal factors alone. In reality, however, one may be able to remove, 
to a large degree, the Lyman continuum opacity due to intergalactic absorbers by making use of a 


tight correlation between Lya and LyC absorption (Inoue & Kamaya 2008). 
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Fig. 4.— Top two panels show the probability distribution function of the apparent mean for 
a given observational sample size iV stack for the high mass (top) and low mass (bottom) sample, 
respectively. The mean is computed by weighting the number of photons produced in each galaxies 
to mimic the stacking of the SED in observations. The true mean of the distribution is denoted in 
each panel. Bottom two panels are the same as the top two panels, but for the subsamples with 
different star formation rates, as indicated in the legend. 
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4. Conclusions 


We have simulated a significant sample of galaxies that are resolved at 4 parsec scales, important 
for capturing the structure of the interstellar medium (e.g., Joung & Mac Low j 2006). We have also 
implemented a much improved supernova feedback method that captures all phases of the Sedov- 
Taylor explosion solution and has been shown to yield the correct final momentum driven by the 
explosion regardless of the numerical resolution (Kimm & Cen 2014). An adequate treatment of 
both these two requirements is imperative, before one can start properly addressing the issue of 
LyC escape, because most of the escape LyC photons escape through “holes” in the interstellar 
medium, instead of them uniformly leaking out in a “translucent” medium. In Kimm & Cen (2014) 
we address the escape fraction for galaxies at the epoch of reionization, to provide the physical basis 
for stellar reionization. 


Here we quantify the distribution of escape fraction for galaxies as a whole, at a range of redshift 
from z = 4 to z = 6. In general, it is found that the LyC escape fraction depends strongly on the 
view angle of the observer and the overall distribution of the escape fraction sampled over many 
sightlines is very broad. The distribution narrows with increasing halo mass or SFR or redshift. 
This broad distribution introduces large sampling uncertainties, when the galaxy sample size is 
limited. For example, a sample of 50 galaxies of halo mass in the range of (10 10 — 10 11 ) M 0 at z = 4 
produces the 2 a fractional range of the escape fraction of 58%-159%. At z = 5 a sample of 20 
galaxies with star formation rate in the range of (0.3 — 10) M 0 yr _1 gives the 2a fractional range of 
the escape fraction is 56%-163%. Our analysis suggests that at least on order of tens of galaxies 
is needed, before one is confident at the 2 a level that the mean escape fraction measured does not 
deviate from the truth by 30-50% at z = 4 — 6 for galaxies hosted by halos of mass in the range 
10 10 - 10 11 M 0 . 
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by the NASA High- End Computing (HEC) Program through the NASA Advanced Supercomputing 
(NAS) Division at Ames Research Center. The research is supported in part by NSF grant AST- 
1108700 and NASA grant NNX12AF91G. 
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(/-2or, /- laif+lcri f+2a)/ (/esc) 


(0.001, 0.004, 3.21, 12.4) 
(0.001, 0.091, 2.34, 6.14) 
(0.091, 0.352, 1.81, 3.48) 
(0.251, 0.509, 1.59, 2.56) 
(0.405, 0.634, 1.43, 2.02) 
(0.580, 0.755, 1.27, 1.59) 
(0.684, 0.823, 1.18, 1.40) 
(0.892, 0.942, 1.05, 1.11) 


(0.001, 0.013, 5.37, 8.79) 
(0.001, 0.003, 3.86, 8.26) 
(0.001, 0.003, 2.58, 7.38) 
(0.001, 0.029, 2.03, 5.68) 
(0.006, 0.151, 1.79, 4.03) 
(0.131, 0.369, 1.58, 2.97) 
(0.277, 0.510, 1.49, 2.41) 
(0.686, 0.820, 1.18, 1.38) 


(0.001, 0.014, 2.89, 9.80) 
(0.001, 0.170, 2.19, 5.18) 
(0.165, 0.443, 1.76, 3.07) 
(0.330, 0.583, 1.54, 2.31) 
(0.476, 0.690, 1.36, 1.84) 
(0.636, 0.797, 1.23, 1.49) 
(0.731, 0.853, 1.16, 1.33) 
(0.907, 0.953, 1.05, 1.10) 


(0.001, 0.001, 2.87, 8.75) 
(0.001, 0.023, 2.21, 5.40) 
(0.034, 0.237, 1.75, 3.70) 
(0.161, 0.395, 1.52, 3.12) 
(0.304, 0.515, 1.38, 2.79) 
(0.468, 0.638, 1.35, 2.10) 
(0.569, 0.711, 1.32, 1.75) 
(0.813, 0.896, 1.10, 1.22) 


(0.001, 0.025, 6.95, 10.1) 
(0.001, 0.009, 5.31, 9.64) 
(0.001, 0.001, 4.33, 9.36) 
(0.001, 0.005, 3.66, 9.15) 
(0.001, 0.025, 2.78, 8.18) 
(0.015, 0.132, 2.31, 6.04) 
(0.063, 0.239, 2.20, 4.65) 
(0.468, 0.692, 1.38, 1.84) 


(/— 2ai f—laif+lcn /+ 2<r)/(/esc) 

10 10 < M vir /M 0 < 10 11 
(0.002, 0.097, 4.57, 14.6) 
(0.011, 0.231, 2.56, 7.45) 
(0.151, 0.470, 1.87, 3.36) 
(0.317, 0.601, 1.59, 2.48) 
(0.470, 0.703, 1.39, 1.92) 
(0.642, 0.805, 1.24, 1.51) 
(0.738, 0.859, 1.16, 1.34) 
(0.910, 0.951, 1.05, 1.10) 

10 9 < Mvir/Af© < 10 10 

(0.001, 0.264, 11.0, 13.4) 
(0.001, 0.155, 10.8, 13.3) 
(0.001, 0.061, 9.42, 13.0) 
(0.001, 0.080, 4.87, 12.3) 
(0.012, 0.185, 2.44, 7.94) 
(0.134, 0.445, 1.78, 3.23) 
(0.320, 0.601, 1.52, 2.29) 
(0.745, 0.866, 1.14, 1.30) 
0.3 < SFR < 10 
(0.002, 0.106, 3.93, 14.3) 
(0.005, 0.272, 2.15, 5.53) 
(0.216, 0.538, 1.62, 2.64) 
(0.410, 0.666, 1.42, 1.99) 
(0.562, 0.760, 1.29, 1.63) 
(0.708, 0.843, 1.17, 1.37) 
(0.787, 0.888, 1.12, 1.25) 
(0.929, 0.964, 1.04, 1.07) 
0.01 < SFR < 0.3 
(0.001, 0.009, 3.74, 12.1) 
(0.001, 0.108, 2.43, 6.76) 
(0.098, 0.388, 1.83, 3.23) 
(0.271, 0.551, 1.58, 2.39) 
(0.434, 0.669, 1.40, 1.91) 
(0.605, 0.781, 1.24, 1.53) 
(0.709, 0.842, 1.17, 1.36) 
(0.900, 0.946, 1.05, 1.11) 
SFR < 0.01 

(0.001, 0.330, 9.99, 12.0) 
(0.001, 0.263, 10.1, 11.8) 
(0.001, 0.097, 9.73, 11.7) 
(0.001, 0.083, 8.51, 11.4) 
(0.011, 0.124, 4.98, 10.5) 
(0.054, 0.261, 2.37, 6.66) 
(0.133, 0.393, 1.96, 4.43) 
(0.585, 0.765, 1.28, 1.64) 


( /— 2a i f- laif+lai f+2u) I (/esc) 


(0.004, 0.162, 3.19, 7.16) 
(0.092, 0.358, 2.41, 4.58) 
(0.311, 0.577, 1.65, 2.63) 
(0.460, 0.682, 1.42, 1.96) 
(0.588, 0.763, 1.28, 1.61) 
(0.717, 0.848, 1.17, 1.36) 
(0.795, 0.889, 1.12, 1.24) 
(0.930, 0.965, 1.03, 1.07) 


(0.001, 0.194, 9.98, 13.1) 
(0.001, 0.132, 9.31, 12.7) 
(0.001, 0.074, 6.97, 12.3) 
(0.001, 0.108, 3.05, 10.4) 
(0.023, 0.231, 1.85, 4.86) 
(0.170, 0.536, 1.49, 2.33) 
(0.419, 0.686, 1.33, 1.77) 
(0.812, 0.900, 1.10, 1.20) 


(0.008, 0.185, 4.74, 9.28) 
(0.062, 0.305, 2.28, 5.49) 
(0.277, 0.555, 1.58, 2.43) 
(0.447, 0.672, 1.39, 1.89) 
(0.583, 0.764, 1.26, 1.58) 
(0.719, 0.846, 1.16, 1.34) 
(0.795, 0.890, 1.11, 1.23) 
(0.931, 0.963, 1.03, 1.07) 


(0.001, 0.025, 3.99, 10.2) 
(0.001, 0.121, 2.02, 5.28) 
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